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The average shape of the Spectral Local Density of States (LDOS) and eigenfunctions (EFs) has 
been studied numerically for a conservative dynamical model (three-orbital Lipkin-Meshkov-Glick 
model) which can exhibit strong chaos in the classical limit. The attention is paid to the comparison 
of the shape of LDOS with that known for random matrix models, as well as to the shape of the 
EFs, for different values of the perturbation strength. The classical counterparts of the LDOS has 
also been studied and found in a remarkable agreement with the quantum calculations. Finally, by 
making use of a generalization of Brillouin-Wigner perturbation expansion, the form of long tails 
of LDOS and EFs is given analytically and confirmed numerically. 
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I. INTRODUCTION 

Recently, growing attention has been payed to the 
structure of the so-called Local Spectral Density of States 
(LDOS) in application to both disordered and dynamical 
systems vifhich exhibit strong chaotic properties (see, for 
example, ) . This quantity, known in nuclear 

physics as the "strength function" , is of special interest 
since it gives information about the "decay" of a specific 
unperturbed state into other states due to interaction. 
In particular, the width of the strength function defines 
the effective "life-time" of the unperturbed basis state. 

Typically, the shape of the LDOS is assumed to be 
of Lorcntzian form (i.e., the "Breit-Wigner shape"), as 
can be analytically derived for sufficiently weak coupling. 
However, in a direct computation of Ce atom |^] it was 
found that at relatively large distances from its center, 
the LDOS has an abrupt decay which is extremely fast 
(even faster than the exponential). This fact, which is 
quite generic, is due to the finite range of the interac- 
tion in the unperturbed energy basis 0||j9). As a result, 
matrix elements of a Hamiltonian describing a realistic 
physical system, decay very fast away from the principal 
diagonal, thus leading to an effective band-like structure. 

Such band structure of Hamiltonian conservative sys- 
tems can be compared to the one known for unitary evo- 
lution operators describing one-dimensional dynamical 
systems under periodic perturbations, like the paradig- 
matic Kicked Rotator Model (KRM) |o|jlll. Another ex- 
ample is an ensemble of Hermitian Band Random Matri- 
ces (BRM), which is used to describe quasi-lD disordered 
models in solid state physics (see, for example, [ p^|Jl3| 
and references therein). The theory of such "standard" 
BRM is now well developed, see review ll^; however, it 



can not be applied, verbatim, to conservative systems like 
isolated atoms, nuclei, atomic clusters, etc. The reason is 
that the Hamiltonians of these latter systems expressed 
in the basis of the reordered unperturbed states, have an 
additional leading diagonal corresponding to the energy 
density of the unperturbed Hamiltonian. Band random 
matrices with such an additional leading diagonal are 
known as Wigncr Band Random Matrices (WBRM) (see 
|,|jl|Jl|,|l|,il3l). Unhke the standard BRM, the theory 
of WBRM is not well developed. On the other hand, 
these matrices are currently under close attention since 
they are believed to provide an adequate description for 
complex systems (atoms, nuclei, clusters, etc), as well as 
for dynamical conservative systems with few degrees of 
freedom, which are chaotic in the classical limit. 

In this paper we consider a specific dynamical model 
of this type, namely, the so-called Lipkin-Meshkov-Glick 
model |18| . In our study we follow the approach devel- 
oped in |W| where the structure of LDOS and eigenfunc- 
tions (EFs) has been numerically investigated in detail 
for the WBRM. The main result of Q which stems from 
direct comparison of LDOS and EFs, is the discovery of 
the so-called " localization in the energy shell" for conser- 
vative systems with chaotic behavior. It is of great inter- 
est to apply the approach suggested in |Q to dynamical 
systems of interacting particles. 

In this connection it may be interesting to remark that 
it is possible to relate specific properties of chaotic eigen- 
states to such observables as the occupation numbers for 
single-particle levels and transition amplitudes (see de- 
tails in [|9lj2^J|l|,||). The above approach |l|,|2^,|l|,|| 
has been developed for the model of two-body random 
interaction, by assuming completely random two-body 
matrix elements. Thus, it is important to extend this 
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approach to dynamical systems of interacting particles 
with a chaotic dynamics. 

The paper has the following structure. In Section II we 
describe the three-orbital Lipkin-Meshkov-Glick model 
and discuss its general properties. The classical limit 
is considered in Section III, where transition to chaos is 
studied in dependence on the strength of the perturba- 
tion. Section IV is devoted to the discussion of general 
properties of eigenstates and spectrum statistics for the 
quantum model. In Section V we numerically investigate 
the structure of LDOS and eigenfunctions for different 
values of model parameters. In section VI, we present 
some analytical and numerical results for long tails of 
LDOS and eigenfunctions. Concluding remarks are given 
in Section VII. 



II. THREE-ORBITAL LMG MODEL 

The three-orbital Lipkin-Meshkov-Glick model [|l8|, or 
in short, the LMG model, is known as some simplification 
of the shell- model of the nucleus. It was introduced also 
to check the validity of approximate many-body tech- 
niques, including the random-phase approximation and 
Bardeen-Cooper-Schrieffer (BCS) theory. The symmet- 
ric states of the LMG model, which we will use in our 
calculations, correspond to collective motions which may 
mimic the collective motion of the nucleus. 

The model has n particles distributed among three 
single-particle orbitals with the same parity and angu- 
lar momentum. Each orbital is fi-fold degenerate. The 
ground, first and second excited orbitals are labeled by 
r — 0,1, 2, and the degenerate states within each orbital 
are labeled by7 = l,2,3,---,ri. The energy of each or- 
bital is denoted by e^. In our calculations, for simplicity 
(without the loss of generality), we will set eo — 0. 

The Hamiltonian of the model is 



H = + \V 



(1) 



where 
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where al^ and fermionic creation and annihilation 

operators obeying the usual anti-commutation relations. 



and parameters A, /^i, ^2, M3, describe the strength of 
the perturbation. 

Hamiltonian (Q) can be expressed in a much simpler 
form. To this end we introduce the two-fermion operators 



Kr 



^ a 



r,s = 0,1,2. 



(3) 



The operators Koq,Kii and ^^22 are number operators 
of the orbitals 0, 1 and 2, and Krs for r 7^ s are par- 
ticle raising and lowering operators, respectively. The 
commutation relations for Krs are 



(4) 



As a result, the Hamiltonian can be written as a function 

of Krs, 



H° = eiA'n + 62^^22 



(5) 



t=i 



where 



yd) = KwKio + KoiKoi 

= K20K20 + K02K02 

= ^21^^20 + K02K12 

V<-^^ =Ki2Kio + KoiK2i. 



(6) 



The nine operators Krs have a very important prop- 
erty, namely, they are invariant under interchange of the 
single-particle-state labels 7. Thus, the Hamiltonian also 
is invariant under interchange of 7 and conserves the per- 
mutation symmetry of the labels 7. This makes possi- 
ble to divide the Hilbert space into subspaces according 
to permutation symmetry. In our quantum calculations, 
we use a subspace composed by the so-called symmetric 
states. A convenient basis \mn > for such subspace can 
be obtained by operating the symmetric raising opera- 
tors Kio and if 20 on the state with all the particles in 
the ground orbital, labeled by |00 >, 



\mn C(m, n)iC™^2"o|00 >, 



(7) 



where C{m,n) is the normalizing coefficient. These 
states |mn > are eigenstates of the number operators Kn 
and K22 with m being the number of particles in the or- 
bital 1 and n the number of particles in the orbital 2. By 
conservation of particles number, = Kqq -f Kn + K22, 
the population of the ground orbital is $7 — m — n. The 
dimension of the symmetric subspace is iV = + l)(ri-|- 
2)/2. Notice that |00 >= 4o ' ' ' ao2«oi |00 • ■ ■ >, and 
therefore, \mn > are antisymmetric under interchange of 
label 7. 

From Eq. (Q) it is seen that the particle raising oper- 
ators KiQ and K20 commute (as a consequence of their 
two-fermionic feature), and therefore, the state \mn > 
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in (|^) is symmetric under the interchange of the order 
of raising particles. In this sense, raising fermions from 
the ground orbital to the two excited orbitals is similar 
to creation of bosons. Indeed, resorting to the generator 
coordinate method approach to the dynamic group rep- 
resentation, a boson representation has been found for 
the symmetric states [ p3|j2^ ]. The above basis states in 
the boson representation are 



mri > = 



{biribiY 



V(m + l)!(n+ 1)! 



100 > 



(8) 



where bj. are creation operators of bosons. The relations 
between the operators Krs and the creation and annihi- 
lation operators of bosons fej and br are 
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KrO^Kl^bl^n-b%-b\b2 
[brMs]^5rs. [br,b,]^[blb\]=Q 



(9) 



for r, s = 1,2. Making use of these relations it is easy 
now to obtain the expressions for the matrix elements of 

^(fi - TO - n)(ri - TO - n - l){m + 1)(to + 2)5m' ,rn+25n' ,n 

< m'n'\K2aK2o\mn >= 

^y {fl - m - n){fl - m - n - l){n + l)(n -|- 2)6m',mSn',n+2 

< m/ n'\K 21 K2o\mn 

^^m{Q - TO - n){n + l){n + 2)6jn',rn-lSn',n+2 

< m'n'\Ki2Kio\mn >= 



y/n{fl - TO - n){m + l){m + 2)S„i>^m+25n 



(10) 



which can also be found from the commutation relations 
&, see H. 

The states \mn > are eigenstates of H° with eigenen- 
ergies 



other states |0j > with the same energy difference d^*^ = 
1^0 - E^\ where = 2ei, d^^) = 2e2, rf^^) = 2e2 - 

ei, d^"*^ = 2ei — £2- The average coupling strength =< 
> can be found by averaging over the non-zero matrix 

elements only. Similarly, for each V^-^^ one can introduce 
(y^l'^y > with the average taken over only the 

non-zero matrix elements of V^-*"^ , respectively. Therefore, 
pit) ^ I flit) ^ natural measure of the strength of 
with respect to the energy distance between the basis 
states coupled by 1/^*^. The parameters ^Xt in Eq. (||) 
are determined by the condition that p*^*^ — 1 for t = 
1,2,3,4, so that, on average, the relative strengths of 
are the same. Under this condition, we have /ii « 
0.0116, ^Ji2 ~ 0.0169, ^Ji^ « 0.0158, « 0.00439, and the 
estimate of the average coupling strength is w w 2.24. 

The global structure of the Hamiltonian matrix is pre- 
sented in Fig. 1, where points represent non-zero off- 
diagonal elements < > of the Hamiltonian H. 
As one can see the matrix is sparse and band-like. More 
precisely, the non-zero elements of the perturbation form 
only eight curves since a basis state > is coupled by V 
to at most eight other basis states. The two inner curves 
result from the contribution of V'^^\ while the two outer 
curves come from V^'^\ Since d'l' = 2.2 w = 2.12, 
the curves corresponding to V'^^'^ and y*^^-* are very close 
to each other and are not separated in the figure. The 
half-band width b depends on Ef and can be analyt- 
ically estimated to be 6 « ?>El/e2 for E^ < E' , and 
b w 2r2-6(£;°-£;')/e2 for E^ > E', where E' « 2^262/3. 
In particular, the maximum width (in the center of the 
band) is bmax ~ 2fl. 

The unperturbed density of states p{E^) of is 
shown in Fig. 2a, and turns out to be in agreement with 
the estimate p{E^) « h/S-'^\ The perturbed density 
of states p{E) is shown in Fig. 2b for A ~ 2.0. For 
a better comparison of p{E) with p{E^), it is conve- 
nient to rescale to the same total energy interval, namely 
p{E) Piy{E) ~ vp{Ev), where v is 



Kin = mei +ne2. 



(11) 



It is convenient to rearrange the eigenstates of in 
order of increasing energy, and we will label them by 



(12) 



Correspondingly, the eigenstates of the total Hamiltonian 
H, also reordered in energy, will be labeled by IV'q > 



H\^a >= Ea\^a > ■ 



(13) 



In our numerical computations on the LMG model, we 
take n = 40; therefore, the dimension of the symmetric 
subspace is A'" = 861. As to the choice of ei and £2, we 
assume, without loss of generality, ei = 1.1 and £2 — 1.61. 

The expressions for in Eqs. (|) and (|l|), show 
that each y^*) couples the basis state \4>i > with only two 



(£^861 — El) 



(14) 



Here, Esei and -Eggi are the highest eigenenergies of H 
and respectively. The ground state energy of the 
perturbed Hamiltonian has also been shifted to coincide 
with the unperturbed one. From the result of such rescal- 
ing, shown in Fig. 2c, one can conclude that the rescaled 
perturbed density of states is similar to the unperturbed 
one. For weaker perturbations, A < 1, the correspon- 
dence is much better. From Fig. 2c one can also see that 
the peak of p{E) is shifted a little towards the center 
of the spectrum. Moreover, for even stronger perturba- 
tions, the peak has been found to be at the center of 
the spectrum. These properties of the global structure of 
the Hamiltonian and of the density of states will be used 
below when we discuss other properties of the model. 
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III. THE CLASSICAL LIMIT 

The classical limit of the symmetric subspace of the 
LMG model can be obtained by two methods. One 
method was used in and consists in the direct study 
of the motion of coherent states in the limit r2 — > oo. 
The other method, which will be used here, is based on 
the boson representation. However, since the boson rep- 
resentation of the symmetric states is obtained via the 
coherent state representation, the two methods are basi- 
cally equivalent. 

In order to obtain the classical limit, we introduce the 
transformation, 

K = \l^ilr - iPr), br = ^^(gr + «Pr), (15) 

for r — 1,2. According to (||) and (p^, qr and Ps obey 
the following commutation rules. 



[qr,Ps 



n 



(16) 



Therefore, the factor plays the role of Planck con- 
stant. Letting the particles number f2 — > oo while keep- 
ing constant the following parameters 



n't = A^t^^^ 



i = l,2,3,4. 



(17) 



one obtains the classical counterpart of the Hamiltonian 



H, 



Hi, 



(18) 



where 



H'd = i{pl + <ll) + i{pl+ll) 

= M'i(gf-P?)(l ~ G/2) +M2('?2 -P2)(l - G/2) (19) 
+ - pDqi + 2g2PiP2] Vl - G/2 



+ - Plh2 + 2gipiP2]v/l-G/2 



2 

-Pl- 



pI 



2{b\bi+blb2)/n < 2. Notice 
depends also on momentum 



with G = qi 

that the perturbation Vd 
variables. 

In order to understand the qualitative properties of 
the classical model, we have plotted the Poincare sur- 
faces of sections at different energies. As in [ p3|j25| ], it 
was found that regular regions of phase space are grad- 
ually destroyed when A increases. However, due to the 



specific form of the classical Hamiltonian Hd in Eq. (19), 
it has been found that the motion on low and high en- 
ergy surfaces can exhibit more chaotic features than on 
the medium energy surfaces. Three typical examples for 



A = 0.9 are shown in Fig. 3. The first figure. Fig. 3(a), 
shows the surface of section at energy E = 10. It can be 
seen that trajectories on this energy surface are chaotic 
except in a small region. The next Fig. 3(b) corresponds 
to the energy E — 39, in the middle of the energy region. 
Here one can distinguish three regular islands. Finally, 
Fig. 3(c) corresponds to the energy i? = 57 in the high en- 
ergy region. Here most trajectories are seen to be chaotic 
with some remaining regular islands. Note that the cen- 
tral region in Fig. 3(c) is energetically inaccessible. When 
A increases to 2.0, it has been found that almost all reg- 
ular islands disappear and the system is nearly totally 
chaotic (Fig. 4). 

The above peculiar behavior is due to the particular 
structure of the perturbation Vd- Indeed, at energy E = 
39, the value of G can be closer to 2 than for the case 
with E = 10. Therefore, due to the terms containing 
(1 — G/2) in (pj|), the perturbation at the energy E — 10 
is stronger than a,t E — 39. For high energies, instead, 
since G is quite close to 2 the derivatives dVd/dpi and 
dVd/dqi can be large and, as a consequence, the motion 
at high energies {E = 57) is more irregular than that at 
the middle ones. 



IV. GENERAL PROPERTIES OF EIGENSTATES 
AND SPECTRUM STATISTICS 

In the above section, we have discussed the classical 
counterpart of the LMG model. In particular, we showed 
that for not very large perturbations , the classical mo- 
tion on low and high energy surfaces is more irregular 
than on the middle ones. In this section, we study some 
general properties of the quantum model, which are re- 
lated to the above classical features. 

In Figs. 5, 6 we show four typical eigenstates \4>a > of 
the total Hamiltonian 77 for A = 0.9 and a = 50 — 53 and 
430-433, respectively, in the basis states \<j)i > (many- 
particle states of the unperturbed Hamiltonian H^). For 
low levels a = 50 — 53, the states \4>a > mainly oc- 
cupy the region i — 200 of the basis states \(f>i >, 
and the expansion coefficients look random in the re- 
gion. For the levels in the middle of the energy spectrum, 
a = 430 — 433, the components of the eigenstates \ipa > 
are mainly distributed in the region i — 200 700, but 
the expansion coefRcients do not appear as completely 
random. For example, the coefRcients < (t>i\^a > for a 
=433 seem to be random, without any structure; instead, 
for a =431, they look sparse and some structure is seen. 
These figures suggest that eigenstates with low energies 
are more chaotic in the region i € [1,200] than those 
with middle energies in the region i e [200,700]. This 
is also confirmed by nearest-level-spacing distributions. 
In Fig. 7 we plot the nearest-level-spacing distributions 
P{s) for eight different regions in the energy spectrum 
of H for A — 0.9. In order to achieve better statistics 
we have diagonalized the Hamiltonian with five different 
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values of A close to A = 0.9 and put together the un- 
folded sequences AEa- As expected, the histograms of 
P{s) for the lowest and highest energy regions are closer 
to the Wigner-Dyson distribution (dashed lines) than to 
the Poisson distribution (dashcd-dot lines), while for a 
in the interval [440,550] P{s) is closer to the Poisson dis- 
tribution. On the other hand, when A increases to 2, the 
level spectrum distribution P{s) become very close to the 
Wigner-Dyson distribution even in the middle energy re- 
gion. 

The above numerical results are related to properties of 
the perturbation V, which are determined by the four op- 
erators KioKio, K20K20, K21K20 and K12K10. So one 
needs to study only non-zero matrix elements of these 
four operators. For example, according to Eq. (10), for a 



fixed basis state \4>i >, there is only one basis state \4>j > 
for which the matrix element < (f)j\KioKio\(j)i > is non- 
zero. Therefore, the non-zero matrix elements of KwKio 
can be regarded as a function of i only. The same is true 
for the other three operators. The dependence of non- 
zero matrix elements of the above four operators on i, is 
presented in Fig. 8. Several features can be seen from this 
figure. Firstly, on average, the non-zero matrix elements 
of ^iKiqKio are relatively large in the low energy region. 
Secondly, apart from the two edges, the average values 
of the non-zero matrix elements of /^2^20-?^20 are simi- 
lar in different energy regions. However, in the middle 
of the energy region, the operator /.t2-f''20^20 has many 
very small non-zero matrix elements. Thirdly, the matrix 
elements of /i3-Hr2i^20 are relatively large, on average, in 
the high energy region. Finally, the variation of the ma- 
trix elements of iiiKi2KiQ in different energy regions is 
not so large as compared to the other three operators. 
As a result, the perturbation is stronger in the low and 
high energy regions than in the middle energy region. 



STRUCTURE OF LDOS AND 
EIGENFUNCTIONS 



quite important for the understanding of generic prop- 
erties of the quantum model. In particular, the LDOS 
shows how the unperturbed state \(j)j > is coupled to 
the exact states \ipa > with the specific energy Ea. The 
width of this function ("spreading width") defines the 
energy range associated with the "life time" of an unper- 
turbed state |0j >. 

The form of the LDOS for band random matrices has 
been analytically studied by Wigner |Q, see also [|. 
Particularly, it was shown that when perturbation is not 
large, the LDOS has the form 



WBw{E~E") = 



r/27r 



{E - k0)2 + r2/4 



(21) 



which is nowadays known as the Breit- Wigner (BW) law. 
Here, F is the half-width of the distribution. For larger 
perturbations, the form of LDOS becomes model depen- 
dent and in the intermediate region can be approximately 
described by a Gaussian distribution p^ . 

Another important quantity is the shape of eigenfunc- 
tions (EFs) 



W^{E") = Y,\C^,fSiE°~Ef) 



(22) 



in the unperturbed energy basis. In our numerical calcu- 
lations of the LDOS and EFs for the LMG model, in order 
to suppress fluctuations, we have taken averages over 200 
of individual distributions in the interval 331 < j < 530 
for the LDOS and 331 < a < 530 for the EFs. The aver- 
aged distributions will be denoted by w{E) and W{E°), 
respectively. Before averaging, we should first express 
Wj{E) and Wa{E'^) with respect to their centroids, re- 
spectively. For the LDOS, the centroid of Wj (E) is just 



Ej — ^ Ea\Ca 



(23) 



In this section we discuss the shape of LDOS and of 
eigenfunctions for the LMG model, we study the classical 
counterpart of the LDOS, and finally we discuss to what 
extent the LMG model can be associated with a band 
random matrix model. 



so that we can express the LDOS as Wj{E — Ej). On 
the other hand, the centroid of Wa{E'^), labeled by Cq, 
is defined by 



J I ^aj I 



(24) 



A. Structure of LDOS and eigenfunctions 

The so-called local spectral density of states (LDOS) 
for an unperturbed state |0j > is defined as 



,(i?) = ^|a,f<5(i?-i?„) 



(20) 



where Ea is the eigenenergy of the perturbed eigenstate 
lipa > and Caj =< 4>j\ipa >■ The function Wj{E), also 
known as the "strength function" or "Green spectra" , is 



and Wa can be expressed as a function of the shift (i?" — 
ea). 

The dependence of the shape of LDOS and EFs on 
the perturbation is presented in Fig. 9. The left column 
gives the LDOS and the right column shows the EFs. 
(Notice that the vertical scale depends on the value of 
A.). The first remark is that the shapes of LDOS and of 
EFs are quite similar when the perturbation is not large 
(A < 0.9). On the other hand, with increasing A, they 
start to deviate from each other. Another result is that 
for not large perturbation, A < 0.9, there seem to be 
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large peaks which are not washed out by the averaging 
process over 200 distributions. In fact, they come from 
dynamical interference (correlation) effects, which will be 
explained in the appendix B. 

The dashed curves in Fig. 9 correspond to the best 
fit to the BW and the dashed-dot curves correspond to 
the best fit to the Gaussian. The fitting was made here 
for the central parts of the LDOS and EFs. Specifically, 
neglecting the long tails, we chose for the fitting only the 
data with w and W larger than 0.01. 

In Fig. 9 it can be seen that for A — 0.5 the central 
parts of both distributions can be well fitted by both the 
BW and Gaussian forms (however, the agreement with 
the Gaussian extends to the region of the tails). For 
stronger perturbation A — 0.9, the distributions can still 
be fitted quite well by the Gaussian form. Finally, when 
the perturbation is very strong, for example A = 2.0, the 
LDOS and EFs deviate, as expected, from both the BW 
and the Gaussian form. 

The difference between the LDOS and the EFs for large 
perturbations A — 2.0, is quite evident. However, one 
should note that the LDOS is plotted in the perturbed 
energy basis while the EE in the unperturbed one. There- 
fore, in order to make the comparison meaningful, one 
should rescale the distribution in a proper way. We used 
the same rescaling as in Fig. 2c. After this rescaling (see 
Fig. 10) they look more similar to each other than in Fig. 
9. 



B. The classical limit of LDOS 

The classical counterpart of LDOS, in short, the clas- 
sical LDOS, labeled by Wci{E — E'^), can be defined as 
the probability that a phase point, which belongs to the 
torus corresponding to the quantum numbers rrij and rij 
of >, has total energy E H. It is expressed as a 
function of the distance {E — Ej) where Ej is the un- 



perturbed energy of the torus 
(p^), in the limit Q ^ oo we have 

2 



According to Eq. 



(25) 



Thus, the torus corresponding to mj and rij is that with 
(Pi -I- ql) = 2mj/n and {pi + qf) = 2nj/VL. 

In analogy to the quantum case, the classical LDOS 
was averaged over 200 different tori. In Fig. 11, we show 
a comparison between the quantum and classical LDOS 
for A = 0.3,0.5,0.9 and 2.0. As discussed in appendix 
B, the LDOS for A — 0.3 depends on strong dynamical 
quantum correlations; as it is seen from the figure, in the 
classical limit these correlations disappear. At the same 
time, the shape of the classical LDOS is, on average, 
close to the quantum one apart from the tails which are 
of quantum origin and are due to tunneling effects. The 
data also show that with increasing the perturbation, the 
agreement becomes better. The main difference for A = 



2.0 is that in the center the LDOS is lower than the 
classical LDOS. This may be related to a result given 
in Ref. |^ that there is a local minimum in the center of 
the LDOS for relatively strong perturbations. 



C. A band random matrix model 

Recently it has been shown that realistic conservative 
systems with chaotic properties can be approximated by 
band random matrices [^]|^,^,^. Therefore, it is very in- 
teresting to check whether the Hamiltonian of the LMG 
model, which does not contain any random matrix ele- 
ment, see Eq. (10), can be associated with an ensemble 
of random matrices. To be as close as possible to the dy- 
namical model, we introduce here band random matrices 
of the form 



Elraii. — Eln 



XVr, 



(26) 



where Hq is the same as for the LMG model, see Eq. (|^) , 
and Vran is obtained by replacing the non-zero matrix el- 
ements of V of the LA/IG model by random numbers with 
Gaussian distribution. The mean value of the matrix el- 
ements (Vran)ki IS zcro and the variance, averaged over 
the non-zero matrix elements, is taken to be the same as 
in the dynamical LMG model < {Vran)ki >=< Vm >■ 

Numerical data for both the LDOS and EFs of iJ^an 
are presented in Fig. 12. Averages have also been taken 
over 200 LDOS and EFs, respectively, in the central re- 
gion of the spectrum, as in the calculations of the LMG 
model. Interestingly, these results are similar to those 
found for Wigner band random matrices that is, for 
small perturbations the central part of the LDOS is of 
the BW form, while in the transition region when the 
perturbation is relatively strong, it can be fitted to the 
Gaussian form; for stronger perturbations it can be fitted 
approximately to the semicircle law 



w{E) 



Rl-E\ 



(27) 



From Fig. 12 and Fig. 9, it can be seen that the shape 
of LDOS and EFs of if™„ for A = 0.3,0.5 and 0.9 are 
much more smooth than the corresponding LDOS and 
EFs of the LMG model. This can be explained by the 
randomness of (K-aTt)fei, see appendix B. Another feature 
is that the central parts of the LDOS and of the EFs 
for A = 0.3 0.9 in Fig. 12 are lower than those for 
the Lipkin model. This also is an effect of interference. 
Finally, comparing the two figures, it can be seen that 
for A < 0.9 the central parts of LDOS and EFs of the 
LMG model are roughly similar to those of Hmn- Thus, 
when perturbation is not very large, the LMG model can 
be associated with the above band random matrix model 
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VI. LONG TAILS OF LDOS AND 
EIGENFUNCTIONS 



In the previous section, we have discussed the central 
parts of LDOS and EFs. In this section we study their 
long tails with the help of perturbation theory since, as 
discussed in appendix A, long tails are always in the per- 
turbative region. 

First we consider the case of small A for which the 
coefficients \Caj \ decrease very fast as \Ea:—Ej \ increases. 
Let us start with the left tails of EFs, « E^, with 



Ea in the middle energy region. From Eq. ( A4 ) , we have 



C1-.|<0,|V..>PH^#^P (28) 



E^-EO 



As indicated in section II, there are only eight basis 
states, that can be coupled with a given basis state > 
by the perturbation V. Denote these states by \cf)i >, 
Z = /i, Z2, • • • ^8 in order of increasing energy. Notice that 
since the energy differences are generally much larger 
than the local level spacings (about 0.1, on average, see 
Fig. 2a), there are many basis states located between 
each two of the above eight states >. Therefore, if 
the term ICafcp decreases fast enough with decreasing k 
in the region E^ << Ea, the component \Cak\ for k = Is 
will be much larger than the sum of the other seven com- 
ponents, then. 



(29) 



Following the procedure given in App. D of Ref. [Q, 
one obtains the following estimate for the left tails of the 
averaged EFs 



(30) 



where ^ = \Ea — -E°|/2e2 and u'^^ is the average strength 
of V^"^^ in the low energy region {V^i^ in Eq. (^9|) is in 
the low energy region). 

For the right tails of EFs, E'j' >> iJ^, similar argu- 



ments lead again to Eq. (30) with u*^^) changed to the 
average strength of V^^^^ in the high energy region, which, 
according to Fig. 8, is approximately equal to the one in 
the low energy region. 

To obtain an expression for the tails of LDOS for |0j > 
with E'^^ in the middle energy region, we assume that in 
the regions of tails the shapes of different eigenstates are 
similar on average, that is 



ri2 



for Ea> -E^^^Ea- E%. Then Eq. (H) gives 



{Ea-Eoy^'^'^ 



(31) 



(32) 



and, for the two tails of LDOS, one obtains the same 
expression as Eq. with £. = \Ea - E°\/2e2 and v^^'^ 
the average strength of F^^-' in the middle energy region 
(vig^ in Eq. ( ^ is in the middle energy region). 

When A is not small, as indicated in appendix B more 
research work is needed in order to obtain an analytical 
expression for the tails of EFs and LDOS. However, we 
can assume that the tails obey a law somewhat similar 
to ( pO[ ) with u^^) changed to v (since when A is not small, 
the tails are determined not only by V''^\ but also by the 
other F'-*-'). Here similar to the small A case, v is the av- 
erage strength of the perturbation in the corresponding 
regions: for the left and right tails of EFs, the average 
should be taken in the low and high energy regions, re- 
spectively, while for the two tails of LDOS the average 
should be taken in the middle energy region. According 
to Fig. 8, the value of v in the low energy region is larger 
than that in the high energy region, so the right tail of 
the EE should drop faster than the left one, while the 
two tails of the LDOS should be similar. 

We have numerically computed the tails of both the 
LDOS and EFs for the case of weak perturbation A = 0.1 
and the results are shown in Fig. 13 in logarithm scale. 
Each point represents an average over 200 states. It can 
be seen that the tails of LDOS and EFs are quite close 
to each other in agreement with the fact that, for small 
A, they obey the same law given by Eq. ( ^0| ) with similar 
values of w*-^''. Also the agreement between the numerical 
results and the analytical prediction is quite good. 

When the perturbation increases to A=0.3, the tails 
begin to deviate from the prediction (^. However, it 
has been found that for A > 0.3 the tails can be fitted 
quite well to the expression ( ^0[ ) with 1; instead of w'-^-' 
and ^ given by 



e = (- 



lOea 



(33) 



where xq is an adjusting parameter. As an example we 
present in Fig. 14 the results for A = 0.9. The LDOS and 
its fitting curve with a:o = 13 are given in Fig. 14a. It 
can be seen that the agreement is quite good not only in 
the long tail regions, but also in the regions quite close to 
the central part. Fig. 14b gives the tails for the EF and 
the fitting curve with xo — 16. The figure also verifies 
the prediction given above that the right tail drops faster 
than the left one for the EF, while the two tails are similar 
for the LDOS. 

For the band random matrix model (p6|), as indicated 
in appendix B, the long tails of EFs and LDOS obey the 
same law given by Eq. ( pO[ ) . This has been confirmed by 
numerical results. As an example, in Fig. 15, we present 
the tails of LDOS for A = 0.9 and the prediction given 
by Eq. (|30| ) . The agreement between numerical data and 
the analytical result is again quite well in the long tail 
regions. 
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VII. CONCLUSIONS AND DISCUSSIONS 

In this paper we have studied the Lipkin-Meshkov- 
Ghck (LMG) model in the many-particle basis of non- 
interacting particles. The main attention has been paid 
to the structure of the local spectral density of states 
(LDOS) in comparison with that of eigenfunctions (EF) 
defined in the unperturbed energy basis. 

Due to its dynamical nature, the properties of the 
LMG model strongly depend on the energy region. 
Namely, for strong enough perturbation A, chaotic prop- 
erties of the model for low and high excitation energies 
are stronger than in the middle of the energy spectrum. 
This fact is explained both on the ground of peculiarities 
of the quantum model and of its classical counterpart. In 
particular, the eigenstates in the middle of the spectrum 
are more regular as compared with the eigenstates for 
low and high energies; correspondingly, the level spac- 
ing distribution is close to the Wigner-Dyson type at the 
edges of the spectrum, unlike in the center of the spec- 
trum where deviations from the Wigner-Dyson have been 
detected. 

One of the main questions addressed in our study is the 
dependence of the shape of LDOS and EFs on the per- 
turbation strength. Numerical analysis has shown that 
for relatively weak perturbations the form of LDOS is 
close to the Breit-Wigner apart from tails. This fact is 
in accordance with several observations for models with 
random interaction described by different random matrix 
ensembles ||l|,^. However, the detailed study of LDOS 
and EFs for the LMG model in the region of not very 
strong interaction, have revealed remarkable correlations 
which can be analytically explained. In fact, these corre- 
lations are due to the dynamical nature of the model and 
they are found to be washed out for stronger interaction. 

With the increase of the perturbation, the form of the 
LDOS changes and for a quite moderate perturbation it 
is well described by a Gaussian. This observation is quite 
interesting in view of recent numerical data for complex 
atoms |I] and heavy nuclei Q|, where the form of the 
LDOS was found to be quite close to a Gaussian. The 
same effect (the change of the form of LDOS from the 
Breit-Wigner to the Gaussian-like) has also been found 
in numerical investigations of Wigner Band Random Ma- 
trices Therefore, our data for the dynamical LMG 
model indicate that the above fact is of quite generic 
nature and occurs also in dynamical models exhibiting 
strong chaos in the classical limit. Finally, when the per- 
turbation is very strong, the LDOS has been found to 
have a quite specific form which is related to the pecu- 
liarity of the model under consideration. 

Another problem is the relation between the shape of 
the LDOS and that of EFs. We have found that after 
a proper rescaling, the shape of the EFs is similar to 
the shape of the LDOS, if the perturbation is not very 
strong. This result indicates the absence of localization 
in the energy shell which has been found in Wigner Band 



Random Matrices 

As was noticed in the shapes of both LDOS and 
EFs have an analogy in the classical limit. In our pa- 
per the relation between LDOS and the corresponding 
classical quantity has been checked for the first time in a 
dynamical model with chaotic classical counterpart (see 
also recent paper |^^). Numerical analysis of the classi- 
cal model has shown that the form of LDOS is close to 
its classical counterpart if the perturbation is not very 
weak. This fact allows to expect that (in semiclassical 
regions) the global structure of both the LDOS and EFs 
can be directly found from the corresponding classical 
model. This is important in view of recent results ^9 22 
revealing the direct connection of the shape of EFs to the 
partition function of isolated systems. 

Of special interest is the question about the applicabil- 
ity of random matrix models to a given dynamical system 
in a chaotic region. There are many examples when full 
random matrices or band random matrices can serve as 
good models for the description of statistical properties of 
spectra and eigenstates of dynamical models. However, 
in these examples the two-body nature of an interaction 
is, typically, not taken into account. In this paper we 
have carefully analyzed the possibility of a random ma- 
trix description of the model in the energy region where 
the corresponding classical system can be treated as a 
chaotic one. Specifically, we have used the same unper- 
turbed part of the Hamiltonian, but with off-diagonal 
matrix elements chosen at random with the same vari- 
ance as in the original dynamical model, keeping zero 
matrix elements which are due to the specific form of the 
interaction. This approach also allows us to reveal to 
what extent the underlying correlations of the dynamical 
model are essential for its statistical description. Nu- 
merical results with such a random model have shown a 
quite good agreement for global properties of the LDOS 
and chaotic EFs. In particular, the shapes of LDOS and 
EFs turn out to remain similar in a large range of the 
perturbation strength. On the other hand, the form of 
LDOS and EFs in the random model turns out to be more 
smooth and do not reveal any regular deviations due to 
quantum dynamical correlations, which, for a weak in- 
teraction, are significant in the dynamical model. 

Finally, the question of long tails of LDOS and EFs 
has been studied in detail, both analytically and numer- 
ically. For this, a generalized approach has been devel- 
oped based on the standard Brillouin- Wigner perturba- 
tion expansion. Namely, the perturbation theory has 
been extended for strong perturbation in the region of 
long tails. This has allowed to find the analytical form 
for the long tails of LDOS and EFs. Numerical data have 
shown a good agreement with the analytical predictions. 
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APPENDIX A: BRILLOUIN-WIGNER 
PERTURBATION EXPANSION OF 
EIGENSTATES 

Here we introduce a generalization of the so-called 
Brillouin-Wigner perturbation expansion ||28| , ^9| , which 
can be shown to be valid even for strong perturbations. 
In particular, we show that long tails are always in the 
perturbative region and this explains also some previous 
resuhs ||,|. 

First, we divide the set of basis vectors, i.e., the eigen- 
states of into two parts, {\(f>i >, i = pi, ■ • ■P2} and 
{\4>i>, « = 1, • • -pi - 1,_P2 + 1, • • -^l- Correspondingly, 
there are two projection operators 



P 



E 

i=pi 



X 



Q = l-P. 



(Al) 



The subspaces related to P and Q will be called in the 
following the P and Q subspaces , respectively. We now 
split an arbitrary eigenstate \ipa > of H into two orthog- 
onal parts 



IV-a >- |i > +1/ >, <t\f>^0 



(A2) 



where \t >= P\4>a > and 
the Schrodinger equation 



>^Q\ipa >■ Applying Q to 
), we readily obtain 



K-HO]|/>-QAl/|Va >, 



then, 



(A3) 



(A4) 



in which the eigenstate \ipa > can be either normalized 
or not. We remark that in the particular case in which 
\t > is chosen as a single basis state > {i ~ a) and 
Ea taken as 



£;„ = s- + <0.|Ay|Va>, 



(A5) 



then the iterative expansion of Eq. (A4) is just the 
Brillouin-Wigner perturbation expansion, in which the 
eigenstate > is not normalized. 



Since our purpose is to study the structure of eigen- 
functions and of LDOS, and not to solve the eigenvalue 
problem, we take Ea in Eq. (A4) as a constant, i.e., its 
exact (unknown) value which can be find out by other 
methods. Then the iterative expansion of Eq. ( A4) gives 



IV'a >=\t> + 



1 



1 



Ea 
rQXV 



HO 
1 



QXV\t > 



Ea-HO^ Ea-H" 

+(1? — T7?^QXVr-'\t> 



QXV\t > 



(A6) 



+( 



EajHO 

E„~H" 



QxvriiJa > 



If the last term on the right hand side of Eq. (A6), de- 
noted by r„, tends to zero when n ^ 00, one has a gen- 
eralization of Brillouin-Wigner perturbation expansion 
(GBWPE), which gives an exact expression for |-0a > 
in terms of It >, Ea, XV and H^: 



+(^^QW|t>+-- 

+ (^-^OAy)^|t>+... 



(A7) 



For example, for a Hamiltonian matrix with band struc- 
ture and perturbation V coupling a basis state to at most 
b other states, if the Q subspace is such chosen that for 
> in the Q subspaces \Ea-E°\ > bVmax, where Vmax 
is the maximum of the absolute value of (AV^j), then T„ 
will approach zero when n goes to infinity. 

For each state \ipa >, there exists a P subspace with 
the minimum number of basis states (denoted by Pmin) 
and correspondingly a Q subspace with the maximal 
number of basis states ( denoted by Qmax) that ensure 
Eq. (A7) to hold. Clearly, from Eq. ([AT]), the components 
of the state |?/;q, > in the Q subspace are in the perturba- 
tive region. Therefore, the expansion of the exact state 
\tl)a > in the basis > is naturally divided into two 
parts: the non-perturbative part \tmin >= Pmin\'4'a > 
and the perturbative part \fmax >= Qmaxl'il^a >■ Tails 
of l^/'a > are always in the perturbative region and can 
be studied with the GBWPE. The same approach can be 
applied to the LDOS, which can also be divided into two 
parts, perturbative and non-perturbative. 



APPENDIX B: A "PATH" APPROACH TO 
PERTURBATIVE EXPANSION 

In order to study the perturbative part of jV'a > in 
GBWPE given in Eq. (A7), we make use of the concept 
of path, in analogy to that in the Feynman's path integral 
theory . First we discuss the case of small perturba- 
tions. In this case, the term \t > can be chosen as the 
basis state > with the smallest value of \Ea — E'^]. 
For an arbitrary > in the Q subspace, Eq. (A7) gives 
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a] =< 4'j\lpa > 



XV, 



- E" Ea~ E" Ea ~ E" 

E 



jk 



kl 



(Bl) 



k,lGQ 



So 



E^ ~ElE^~ E° 



We now term a sequence j — > ki 



kn-l 



a path of q paces from j to i, if the matrix elements of 
the perturbation Vkk' corresponding to each pace is non- 
zero. Attributing a factor XVkk' /{Ea — E'l) to each pace 
k k' , the contribution of a path to Caj is the product 
of the factors of all its paces. Then , the g-th term of Caj 
on the right hand side of Eq. (|B1[) can be rewritten as 



(B2) 



where s indicates the paths with q paces from j to i and 
fq,s(j — > i) is the contribution of the path s to Caj- 
Defining A{i i) as 



(B3) 



9=«o 



where qg is the number of the paces of the shortest path 
from j to j, we have 



Ca 



(B4) 



LMG model it is especially useful since the perturbation 
V couples a basis state to at most eight others as shown 
in Fig. 1. In fact, we have a simple method to find out all 
possible paths from j to i. First we should write the basis 
vectors |0i > in their equivalent forms \mi,ni > with rrii 
being the number of particles in the orbital 1 and rii the 
number of particles in the orbital 2. Let us denote the 
changes of m and n along a path by Am = — uij and 
rij, and the changes of m and n for a pace 
kl to k' — by drrir — mk' — ruk and 

Uk- Then, if the path has q paces, we can 



Hi 



An 

r from k = 
Srir — rik' — 
write 



Am = Sr 

r=l 

q 

An = (5n 



(B7) 



Since the interchange of r does not influence the sum, 
Eq. (B7) is also satisfied by some other related paths. 
According to Eq. (p^), there are only eight possible pairs 
of 5m and 5n 



6m = ±2, Sn = for V^^^ 

5m = 0, 5n^ ±2 for 
(5m = Tl, 5n = ±2 for V^^) 
5m = ±2, 5n = =f1 for V'~^\ 



(B8) 



for small W . 

For the general case, the term \t > cannot be chosen 
as a single basis vector, but must be a superposition of 
some basis vectors. Writing 



P2 

E 



ti\(i}i >, 



Caj —< 'Pjli^a > can be expressed as 



P2 



(B5) 



(B6) 



For Hamiltonians with band structure, paths have 
qualitatively different features for small and large pertur- 
bations. When the perturbation is small, the difference 
(p2 — Pi), see Eq. (|AlD, of the Pmin subspace is smaller 
than the band width h. A path which start from j < pi, 
can reach points k larger than p2, in other words, a path 
can "cross" the P region. Thus, the difference Ea — E^ 
can be both positive and negative. This is important for 
the parameters we have chosen in section II, which make 
all the non-zero matrix elements of V positive. When A 
is so large that p2 — pi is larger than b, paths starting 
from j < pi can no longer reach k > p2, i.e., can not 
"cross" the P region. Hence the values of Ea 
either always positive or negative for a path. 

The above concept of path is very useful in studying 
Hamiltonians with band and sparse structures. For the 



E^ are 



therefore, in order to find out all possible q-pace paths 
from j to i, one should just find out all possible combi- 
nations of q pairs of the possible 5m and 5n in Eq. (BS) 
that satisfy Eq. (B7) (if all the intermediate points are 
in the Q subspace) . 

When A is small, the term |i > in Eq. (AT) can be 



chosen as one basis vector >, and the expanding co- 
efficient Caj =< (pjli^a > is A{j i) = J2q Eq, sce Eq. 



( B4 ) . As an example, let us consider the Caj for which j 
is determined by Am — mi — mj — 2, An — rii — rij — 0, 
i.e., the two unperturbed states > and \4>j> are cou- 
pled directly by the term KiqKiq in Eq. (o). Denote 
the number of the q-pace paths with positive contribu- 
tions fq^s in Eq (see Eq. (B2)) by iVq+ and the number 



of paths with negative /g,s by 7Vg_. For the shortest 
path with one pace, A^i+ = IjA^i- — 0. There is no 
path with 2 paces. For q — 3, there are 12 paths (e.g., 
K02K02 K20K20 KioKio), N3+ = 11, TVs- = 1, 
and for g = 4, we have A'4+ = 11, = 7. Clearly, the 
number of short paths with positive fq,s is much larger 
than that of short paths with negative fq^s- Therefore, 
I Caj I is quite large and the dynamical interference effect 
is strong. Similarly, for the other seven \(j>j > that can be 
coupled directly to > by the other terms of V, \Caj\ 
are also quite large. On the other hand, for \(j>j > which 
are not coupled directly to > hy V, \Caj\ are smaller. 
In fact, the two peaks beside the main peak in fig. 9 for 
A = 0.3 come from Caj for the eight states \(j)j > coupled 
directly by V to the state >. 
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For the band random matrix model discussed in sec- 
tion V, since the non-zero matrix elements of Vran have 
random signs, we have Nq+ ~ Nq^, and the dynamical 
interference effect should be much smaller than that in 
the LMG model. 

Finally, let us give some discussion about long tails of 
LDQS and FFs when X is not smpill Tn this 



"asp, wp 



[2] 
[3] 
[4] 
J5I 



shoilld use F.qs 



and (Rfi) Tn the regimi of long 



tails, for each A{j — > i) the main contribution comes 
from the term Fg„. To see this, we give an estimate of 
Fg{j —^ i). Since one basevector can be coupled to at 
most eight others by V, the number of possible paths 
with q paces from j to i is Af? where < M < 8. Then 
for long tails we have 



Mv 



(B9) 



where v is the average coupling strength of the pertur- 
bation V and for large g, E^ is approximately equal to 
E^. When \Ea — E^\ is large enough, the value Fq{j i) 
decreases very fast with increasing q and the main term 
is the one with the shortest path (go paces) from j to i. 



Therefore, Caj can be estimated as 



(BIO) 



For the LMG model, due to the dynamical interference 
(correlation) effects as discussed above, Fqg{j i) for 
different i may be equally important. In fact, for large 
go, although the number of paces is large, the number of 
paths is also large and some of the paces at k may have 
relatively small values of \Ea — E^\. Thus, it is difficult 
t o ob tain an analytical expression for the tails from Eq. 
(|l^. 

For the band random matrix model ([2^), the signs of 
the non-zero matrix elements of Vran are random. There- 
fore, as discussed above, dynamical interference effects 
are not so large as in the LMG model, and the largest 
term on the right hand side of Eq. (BIO) is the one with 
the smallest number of paces go- Since dP'^ is the largest 
among the four d^*^ , the shortest path from j to i mainly 
consists of paces resulting from F^^). Then Eq. (|9|) also 
holds and the same expressions for the long tails of EFs 
and LDOS as in (BOh can be found. 
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FIG. 1. Global structure of the Hamiltonian matrix. 
The points represent the non-zero off-diagonal elements 
< > of the Hamiltonian H . 
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FIG. 2. (a) Unperturbed density of states of H° . 

(b) Density of states p{E) of the system H for A = 2. (c) 
The rescaled density of states pv{E) for \ = 2 (solid his- 
togram) compared with the unperturbed density of states 
piE") (dashed histogram). 



FIG. 3. Poincarc surface of sections on the {q2,P2) plane 
with gi = for A = 0.9 and (a) E = 10, (b) E = 39 and (c) 
E = 57. 



FIG. 4. Same as in Fig.3 for A = 2 and (a) E=A, (b) E=A3 
and (c) E=62. 



FIG. 5. Four typical eigenstates \tpa > of the Hamiltonian 
H ioT X = 0.9 in the basis \4>i >; (a) a =50, (b) a =51, (c) a 
=52, (d) a =53. 



FIG. 6. Same as in Fig.5 for (a) a = 430, (b) a = 431, (c) 
a = 432, (d) a = 433. 

FIG. 7. Histograms of the nearcst-lovcl-spacing distribu- 
tions P{s) for eight energy regions of the Hamiltonian H with 
A = 0.9. (1) 1< a <110, (2) 110< a <220, (3) 220< a <330, 
(4) 330< a <440, (5) 440< a <550, (6) 550< a <660, (7) 
660< a <770, (8) 770< a <861. The dashed and dashcd-dot 
curves represent the Wigner-Dyson and Poisson distributions, 
respectively. Each histogram was obtained by diagonalizing 
five different Hamiltonians with values of A close to 0.9. 

FIG. 8. Non-zero matrix elements of the four operators (a) 
l^iKioKio, (b) 112K20K20, (c) /i3Jf2iJir2o and (d) p,iKx2Kw. 
Each operator couples a basis state 10, > to at most only one 
other basis state \4>j >. 

FIG. 9. Averaged LDOS (left column histograms) and 
eigenfunctions (EF) (right column histograms) for the LMG 
model with different values of A. Dashed and dashed-dot 
curves are fitting curves to the BW and Gaussian forms, re- 
spectively. 

FIG. 10. A comparison between the rescaled LDOS Wsc{E) 
(sohd histogram) and the EF W{E°) for A = 2 (dashed his- 
togram) . 

FIG. 11. Classical counterparts of the LDOS (dashed-dot 
curves) for (a) A = 0.3, (b) A = 0.5, (c) A = 0.9 and (d) 
A = 2. For comparison the corresponding LDOS are also 
plotted (solid histograms). 



FIG. 12. Same as in Fig.9 but for the band random matrix 
model Hran- The solid curve for the LDOS at A = 2.0 is the 
fitting curve to the semicircle law Eq. (27) with Ro w 23.9. 

FIG. 13. Numerically computed LDOS (circles) and EF 
(triangles) for A = 0.1. Each value represents an average 
over 200 states. The dashcd-dot curve gives the analytical 
prediction (30) (with the appropriate normahzation factor) 
for the long tails. 

FIG. 14. (a) LDOS (circles) for A = 0.9 and the fitting 
curve (dashed-dot curve) given by Eqs. (30) and (33) with 
xo = —13. (b) Similar to (a) for the EF (circles) with 
xo = —16. 

FIG. 15. LDOS (circles) of Hran in Eq. (26) for A = 0.9 
and the theoretical prediction (30) (dashed-dot curve) for its 
long tails. 
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